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' Using tlie variational formula for operator product coefficients a method 

for perturbative calculation of the short- distance expansion of the Spin- 
psj . Spin correlation function in the two dimensional Ising model is presented. 

^ . Results of explicit calculation up to third order agree with known results 

Q>^ I from the scaling limit of the lattice calculation. 



Oh! 



X 



It has been nearly fifty years since exact solvability and existence of continuous 
phase transition in two dimensional Ising model was shown.[] Ever since then, thanks 



H ' to equivalence of statistical mechanics models and quantum euclidean field theories 

with an ultraviolet (UV) fixed point, the Ising model has been an excellent theoretical 
playground for addressing many field theoretical issues, such as the existence and 
proper definition of continuum field theories @ and the coordinate space description 
of the renormalization group (RG). In particular, the existence of a second order phase 
transition in the Ising model is equivalent to the existence of a UV fixed point for the 
corresponding quantum field theory of a massive Majorana fermion. Understanding 
the structure of the short-distance expansion of correlation functions of the lattice 
model in the scaling limit by field theoretical methods amounts to renormalization 
group analysis and perturbation theory near the UV fixed point 0. 



*Supported in part by D.O.E. contract DE-FG03-91ER 40662 Task C. 

^For the Onsager solution of the Ising model, see ref. and references therein. 
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In the spirit of classic works of Wilson on the renormalization group and operator 
product expansion (OPE) 0, a new framework for the study of the theory space of eu- 
clidean quantum field theories has been introduced[0. This formalism is particularly 
useful for elucidating the structure of the short-distance expansions of the ra-point 
functions of a renormalizable quantum field theory around a nontrivial fixed point. 
As a first application of this formalism a complete renormalization group analysis of 
the short- distance expansion of spin-spin correlation function in the scaling limit of 
the Ising model has been given[^]. The RG analysis determines the short-distance ex- 
pansion of the spin-spin correlation function only up to RG invariant functions. The 
calculation of these functions requires a perturbative calculation based on a concrete 
regularization scheme. The purpose of the present paper is to show how this calcu- 
lation is done. To this end we show how the recursive use of the variational formulas 
for the correlation functions and operator product coefficients |0] leads to a man- 
ifestly finite calculational scheme for the short-distance expansion of the spin-spin 
correlation function in the Ising model. 

The scaling limit of the two dimensional massive Ising model is characterized by 
the theory of a free massive Majorana fermion which is parametrized by the cosmolog- 
ical constant and the mass m. The cosmological constant only provides an additional 
constant free energy density, and therefore none of the correlation functions of local 
operators depend on it. Therefore the theory is parametrized by one real parameter 
m. The mass has scale dimension one and satisfies the canonical RG equation. Note 
that mr is RG invariant. The operator conjugate to m, denoted by Om, has the 
following expectation value determined by RG analysis |Q 

{Om)m = -Pim\n{mc), (1) 

where /3i is the anomalous dimension of Om and c is a constant of integration. At 
the fixed point this operator corresponds to the energy density operator S = ^{itpip) 
where ip is a free massive Majorana fermion. Our particular choice of normalization 
fixes /?i = —77-. 

The main tool in this scheme is the variational formula which realizes the deriva- 
tive of correlation functions with respect to m by the volume integral of the same 
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correlation functions with an Om insertion. The necessary subtractions are made by 
introducing the operator product expansion from the outset ||]. By counting scaling 
dimensions the only necessary subtraction in the present calculation is 

Jl>r>e 

= ;^ln(e) +c„„^, (2) 
/vr 

where C'^m"'" is the singular part of the coefficient of the OPE of two energy density 
operators and Cmm^ is the corresponding finite counterterm. Using OPE and the 
variational formula for the correlation functions a manifestly finite variational formula 
for OPE coefficients can be derived^. We show how this VF leads to a manifestly 
finite scheme for the perturbative calculation of the corrections to the scaling behavior 
of the spin-spin correlation function. 

Consider the short distance expansion for spin-spin correlation 

(a(r)a(O))„ = a/(r;m)(0,)^, (3) 

where the sum over b is sufficient to include scalar composite operators, which are 
not total derivatives. Because of the operator algebra of the critical Ising model, 
06 G {[1], To calculate the corrections to scaling behavior of the OPE coefficients 
to third orderQ in mass it is sufficient to truncate the sum over b toQ 

(a(r)a(O))^ = C^^\r;m) + C^r{r;m){Om)m. (4) 

The zero-th order is just the well-known scaling limit. In first order the following 
finite and truncated first order VF for the OPE coefficients can be derived]^, § 



R'^OO I jR'>r' 

-C^r{r) lim ( / rfV (f(r')^(0)) + An^'(e))j. (5) 

^Up to scale dimension 3 the only non-derivative scalar operators are 1 and f j|, p^ . 
^For a general proof of truncation of the variational formula to arbitrary order see ref. M. 
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The free massive Majorana fermion representation for the critical Ising model |T0 



makes the calculation of the energy-density correlation functions trivial. The ex- 
pression of cr(r) in terms of the fermion field, on the other hand, is non-local. This 
makes the calculation of the correlation functions, which involve the spin operator, 
non-trivial^. In two dimensions, all correlation functions of conformally invariant 
field theories are calculable in principle Moreover, in the case of the critical Ising 
model there are additional structures which give a systematic way of recursively cal- 
culating the correlation functions [[1^] needed in the main terms of the VFs {i.e., the 
terms denoted by Ais in what follows). In the following we use the results, which. 



using the techniques of |T^, can be easily derivedg. For example, in the case of VP 
(I) we have 



where 

Mi{U') = [ r, (7) 

^ ' Ju'>u' 47i\u'\\u' - ly ^ ' 

Si£{U')=\mv ( I f,""' +— ln(e)). (8) 

The repeated use of variational formula leads to expressions for higher order 
derivatives of Caa^ evaluated at m = 0. In deriving these expressions up to third 
order we should keep the following two points in mind. First, since we only con- 
sider VP for correlation functions which involve energy- density and spin operators, 
the knowledge of the subtraction (|^) is sufficient. Second, in second and third order 
VFs, where the arguments of two inserted energy-density operators, one coming from 
a lower order VP, approach one another, we do not need any subtractions at all. 
This is due to the connectedness of the insertions, i.e., we have to consider partially 
connected correlation functions with respect to all insertions. Therefore, by taking 
the derivative of VF and using the first order VP for correlation functions, we 
could derive the following second order finite and truncated VF 

■*This is the main reason that, even though the critical Ising model has a representation in terms 
of a free massless Majonara fermion, it is still a non-trivial fixed point theory. 
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^2 a.'(r) = lim lim ( / rfV / dV {Sir")S{r')(7ir)aU])Y 

" R'^OO R"^00 17h'>,' jR">r" 



lim ( / rfV {S{r')Sm + A„^^(e)) + (i?' ^ i?")! j, (9) 



The effect of the terms corresponding to d'^{£) vanish because of the specific order of 
the hmits and the OPE algebra of the critical Ising model. We emphasize that this 
equation is finite^. Using the explicit expression for the necessary four point function 
we have 



d'LCa„'^{r) = lim hm 



,(10) 



where 



u'>u' ATr\u'\\u' — 1| Ju">u" A'k\u"\\u" — 1| 

m' + u" - 2mV'P - A\u'\ |m' - 1 1 \u"\\u" - 11 



In complete analogy to second order we can derive higher order finite and truncated 
variational formulasH]. All the necessary critical correlation functions can be calcu- 



lated with the aid of the recursive relations of ref.[|12|. As a result we have two types 
of integrals to handle. The calculation of the subtraction integrals iS's are straight- 
forward, due to the rationality of the integrands. But, the main integrals A^'s are in 
general quite hard and have been one of the main obstacles for the field theoretical 
calculation of the short-distance expansion of the spin-spin correlation function. Here 
we give a general method for calculating these integrals]^. 

Taking the IR limits to infinity has a very important calculational advantage; 
usually it is much easier to calculate integrals for infinitely large cutoffs rather than 
finite cutoffs. Taking advantage of this fact the integral M.i could be calculated 



^For more detailed explanation of this subtle point see ref. 
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and involves elliptic functions. But, the higher order integrals are very difficult to 
calculate directly. The main reason for this difficulty is the appearance of one copy 
of the integrand of A^i, for each integration. Once we analyze the structure of A^i 
carefully, the analysis of the higher order integrals follows. 

The separation of the two poles of the integrand of A^i makes the integral asym- 
metric respect to the IR and the UV behavior; we do not have any UV divergences 
in the main term. Therefore we have to look for a transformation which symmetrizes 
the integral with respect to the IR and the UV limits. We now demonstrate that the 
following conformal transformation does the job 




hence, the name symmetrization transformation. Under this transformation, for large 
cutoff, the A^i becomes 

After applying the symmetrization transformation, the integral in is trivial 
and for Cmm"*" = - In 2 we get the following results 



a."^(r) = Trri, (14) 
dmCaa^{r) = r"3lnr. (15) 

The logaritmic dependence on r in ( P3[ ) could be also derived from solving the RG 
equation p. In addition, using the VF for {Om)m, and the representation of Ising 
model in terms the theory of the free massive fermions it is possible to calculate the 
complete relation between c in (|I|) and Cmm"*" The final answer for our choice of 
Cmm"*" is c = ^5 where 7 is the Euler constant. This completes the specification of our 
choice of the conjugate operator. The ffist non-trivial information will come out of 
the second order calculation. 

The calculation of involves non-trivial integration of elliptic functions and is 
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extremely liard.Q. But using the symmetrization transformation, in the limit of large 
cutoffs, this integral is calculable with standard methods and we get 



lim hm M2{U',U") = l. (16) 



Therefore, the main term in (|T0|) does not have any IR divergences, and we have 



dmC.rir) = 0, (17) 

7 

dlC^.'ir) = Z. (18) 



Note that, (|17|) is a result independent of the choice of any finite counterterms. In 
ref.[^ this was a result expected from comparison to the exact result of ref. 0. 

The third order calculation also relies on the symmetrization transformation. Just 
like the second order calculation, the result of applying this transformation is an 
integral which can be handled by standard techniques. We finally get 
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dlC^.-^ir) = (19) 

Ca.'(r) = --Inr. (20) 
lb 

Contrary to (|T^, the vanishing result for the constant part of (^) is due to the 
specific choice of the finite counterterm Cmm"*"- We can now calculate the short-distance 
expansion on the spin-spin correlation function to third order in the mass. In complete 
agreement with the scaling limit of the exact lattice results of ref . P] we have 



(a(r)a(0))„ = 4(l + hnitleyS) + ^t' + l-tHn{\t\ey 8) + 0{tHn' \t\)), (21) 
ri ^ Z lb o2 ^ 

^Dotsenko applied his conformal field techniques to calculate this integral by clever contour 
deformations and analytic continuation methods. |^ But, as he also points out, generalizing his 
technique to higher orders is very difficult. 
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where t = mr. This is the first systematic field-theoretic calculation near a non-trivial 
fixed point. 

Note that our method is general, and in principle, there seems to be no difficulty in 
continuing the calculation to higher orders. A complete proof of the truncation, as well 
as a systematic study of the structure of the higher order calculations, applicability of 
the symmetrization transformation to arbitrary order, and details of the calculations, 
is presented elsewhere H]. We believe that our methods should also be applicable for 
perturbative calculations around other non-trivial fixed points corresponding to more 
complicated minimal models classified in ref. |T^. It is important to note that the 
proof of the finiteness and the possibility of truncation is applicable to all models of 
this type|^. Clearly, there is much important work to be done in this direction. 

We would like to thank H. Sonoda for suggesting the problem, his valuable com- 
ments, and constant encouragement. 
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